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ABSTRACT 

A number of algorithms for recursive computation of the 
discrete cosine transform (DCT) have been developed 
recently. This paper presents a new method for computing 
the discrete cosine transform and its inverse using 
Bakhvalov’s algorithm, a method developed for evaluation 
of a polynomial at a point In this paper, we will focus on 
both the application of the algorithm to the computation of 
the DCT -I and its complexity. In addition, Bakhvalov’s 
algorithm is compared with Clenshaw’s algorithm for the 
computation of the DCT. 

1. INTRODUCTION 

The discrete cosine transform (DCT) is axnmonly used in 
signal processing, signal compression, image processing, 
and communication systems [1]. Researchers have 
developed efficient algorithms for computing the forward 
and inverse discrete transforms. These algorithms are 
similar to the fast Fourier transform (FFT) algorithms for 
computing the discrete Fourier transform (DFT). 
Specifically, they are more appropriate for minimizing the 
complexity of computations [2, 3], 

Wang and Hunt have shown that for Markovian signals 
with correlation coefficient of less than 0.8, the symmetrical 
DCT (DCT-I) has better performance than the DCT -II [4]. 
In addition, the DCT-I has been used to compute the real 
discrete Fourier transform [5]. The essential computation in 
all these transforms is of the form 

N 

Y ( k ) = ^ y (w) cos (Jem / AO (1) 

for k = 0,1,. ..,W. 

This paper will focus on the recursive computation of 
the symmetric discrete cosine transform (DCT-I) using a 
computational method developed by Bakhvalov for the 
evaluation of an Nth order polynomial at a point [6]. This 


method developed for polynomial evaluation is more 
accurate numerically than other methods in certain cases 
[6]. The back-to-front computational method computes the 
sum of a sequence of numbers by adding small numbers 
first to minimize the round-off errors. Kronsjo [6] presents 
an upper bound for the round-off error. Bakhvalov’s 
method requires additional adders and dividers by 2 
compared with Clenshaw’s methods [5, 6]. 

Since both Bakhvalov’s and Clenshaw’s algorithms are 
recursive, they are easily implemented on parallel 
computing platforms as well as VLSI and Held 
Programmable Gate Arrays (FPGAs). 

2, Computation of the Discrete Cosine Transform 

The symmetric discrete cosine transform (DCT-I) [1, 7] of 
N+l data points y (0), y (1),. y (N) is given by 

N 

Y(k) = (2/ Nf n a k y^a„y(n) cosQcrm / N) (2) 

«=0 

and its inverse is given by 

N 

y(n)={2/Nf n a„ y £a k Y(k)cos(?ikn:/N) (3) 

fork,n = 0, 1, N, 
where 



when p = 0 or N and 

a p =l 

when 

p&Oor p . 


3. BAKHVALOV’S ALGORITHM FOR EVALUATING 
POLYNOMIALS 




Figure 3: Recursive structure of algorithms 

4. (COMPUTATION OF OTHER FORMS OF DOTS 

The most commonly used DCT is the DCT-EL Wang eL aL 
[8] showed that the DCT of N data points y( 0), y(l), 
y(N-l) can be expressed as 

N-i 

¥(k) = (2 /N) y k y £ i y(N-l-n)cosMn+if2)n/ N)] (7) 

n - 0 

N~l 

y(n) = (— l) n ^r(iV-l-t)sin[(fc+lXn+l/2)7r/Aa (8) 
k=0 

for k,n =0,1 N-l, 

where y k ~ Vi if k=0 and 1 otherwise. 

Using the relationships of equations (9) and (10), 
equations (7) and (8) can be put in the form of equation (2), 
the symmetric DCT form [2,9]. 


constants, ramps, parabolas, exponentials and random sets 
(uniform distributions). Trials were run on each set using 
256, 512 and 1024 data points and several different ranges 
for the amplitudes. 

We used the sorted direct implementation of 
equation (3) as the reference for numerical accuracy. The 
product terms in the direct implementation of the DCT were 
sorted and added from smallest to largest values to 
minimize round-off error when computing the partial sums. 

Table 1 shows the inputs with N data points, 
y(0), y(l), y(2) . . .. Y(N-1) for which the numerical accuracy 
of Bakhvalov’s and Clenshaw’s methods were computed 
using Matlab, where u(n) is the unit step sequence. The 
entries under Bakhvalov indicate the cases where this 
method was more accurate than Clenshaw’s. Similarly, the 
entries under Clenshaw indicate the cases where this 
method was more accurate than Bakhvalov’s. Both 
methods were within approximately one percent of each 
other. It is difficult to know in advance which method will 
produce more accurate results. 

Both Bakhvalov’s and Clenshaw’s methods 
proved to be quite accurate. One method occasionally 
proved to be superior for a given group of data sets, for 
example Clenshaw’s method when used on a ramp and 
Bakhvalov’s for unit steps of length 1024. 

In addition, we tested 50 random input 
sequences of length 256, 512, and 1024 using uniformly 
distributed (0,1) values. The MSE for Clenshaw’s was 
lower in 27 cases for length 256 and 32 cases for length 
512. The MSE was lower for Bakhvalov’s method in 29 
cases for length 1024. 

6. SUMMARY 


(N+l)/2 

^ y(n) cos w(k + 1 / 2) = 

*=i 


(V- 1) / 2 

cos(>v/2) ^ y(n)cos(wk) 
*=o 


(tf+l)/2 

y(n) sin w(k + 1/2) = 


*=i 


(N+ 1)/2 

sin( w / 2) ^ y (n) cos( wn ) 
*= o 


We presented a new method for the computation of the 
recursive symmetric discrete cosine transform based on 
Bakhvalov’s method for evaluating a polynomial at a point 
It is a back-to-front computational method and easily 
implemented on parallel platforms and hardware. When 
) compared with Clenshaw’s method., the recursive part of 

both algorithms has the same complexity. Both 
Bakhvalov’s and Clenshaw’s methods can be used to 
compute other forms of discrete cosine transforms. In 
addition, we compared the numerical accuracy of the two 
(10) methods. 
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